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Abstract 

The electronic and magnetic properties of spin-1/2 pyrochlores R2V2O7 were investigated on 
the basis of density-functional calculations. Contrary to the common belief, the spin-1/2 V i+ ions 
are found to have a substantial easy-axis single-ion anisotropy. The \D/J\ ratio deduced from 
the magnon quantum Hall effect of LU2V2O7, where J is the nearest-neighbor spin exchange and 
D is the Dzyaloshinskii-Moriya parameter, is much greater than the value estimated from our 
calculations (i.e., 0.32 vs. 0.05). We show that this discrepancy is due to the neglect of the single- 
ion anisotropy of the V 4+ ions, and the negative magnetoresistance observed for R2V2O7 arises 
from a new mechanism. 

PACS numbers: 75.30. Gw,71.15.Rf,71.20.-b,75.10.-b 
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During the past two decades the magnetic properties of pyrochlore oxides of the type 
Rj + M2 + 07 (R = rare-earth; M = transition metal) have been extensively studied^, due 
largely to the spin frustration in the pyrochlore lattice that results when the nearest-neighbor 
(NN) spin exchange is antiferro magnetic (AFM) 1 . The vanadate pyrochlores R2V2O7 (R = 
Lu, Yb, Tm, Y) are unique because they are ferromagnetic insulators 3 - - - contrary to the 
observation that ferromagnetism leads usually to metallic character. Furthermore, LU2V2O7 
exhibits a negative magnetoresistance (MR) as high as 50% just above the Curie temperature 
Tc = 73K under magnetic field of 5 T— . Recently, LU2V2O7 is found to exhibit a magnon Hall 
effect (i.e., the anomalous thermal Hall effect caused by spin excitations) 1 ^. The explanations 
presented for these observations raise fundamental questions. Namely, the MR of LU2V2O7 
was suggested to be caused by polaron mediation as found in Tl 2 Mn 2 07 9 . However, this 
possibility seems remote, because the diffuse 6s orbital of the Tl 3+ ion is believed to assist 
the polaron formation in Tl 2 Mn 2 07 whereas the Lu 3+ ion in LU2V2O7 does not have such an 
extended s orbital. In the magnon Hall effect 1 ^, the Dzyaloshinskii-Moriya (DM) interaction 
is considered to play the role of the vector potential as in the intrinsic anomalous Hall effect 
in metallic ferromagnets. In the fitting analysis of their experimental data, Onose et al. 
found it necessary to use the \D/J\ ratio of 0.32. This ratio is unusually large because the 
DM interaction is a consequence of spin-orbit coupling (SOC) so that the typical \D/J\ ratio 
is expected to be smaller than 0.1. Therefore, it is important to quantify the magnitude 
of the DM term. In this Letter, we probe these questions by studying the electronic and 
magnetic properties of R2V2O7 on the basis of density functional calculations to find that the 
spin- 1/2 V 4+ ions have a substantial easy-axis single-ion anisotropy contrary to the common 
belief, the neglect of this anisotropy can lead to an unusually large \D / J\ ratio, and R2V2O7 
exhibits a new type of negative MR mechanism. 

Our calculations are based on the density functional theory (DFT) plus the on-site re- 
pulsion (U) method^ (DFT+U) within the generalized gradient approximation 1 ^ on the 
basis of the projector augmented wave method 14 encoded in the Vienna ab initio simula- 
tion package 1 ^. The plane-wave cutoff energy was set to 400 eV. Careful convergence test 
were performed and the total energy was converged to 10~ 6 eV. In the following, we report 
results obtained with U = 3 eV and J = 1 eV on V, but the dependence of our results on 
U and J will be also discussed. It was found 16 that Y2V2O7 has magnetic properties very 
similar to those of LU2V2O7. Thus, our calculations focused mainly on Y2V2O7 using its 
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experimental structure 1 ^. The structural optimization was found to have a negligible effect 
on its magnetic properties. 

In the V2O7 framework of Y2V2O7 (space group Fd3m), the VOq octahedron containing 
V 4+ (3d 1 , S = 1/2) ions share their corners such that the V 4+ ions form corner-sharing V4 
tetrahedrons [see Fig. [0(a)]. below T c ~ 70 K. The V0 6 octahedron of R 2 V 2 7 are slightly 
distorted (i.e., axially-compressed slightly) from the ideal shape. Under the trigonal crystal 
field, the £23 states are split into the lowest a\ g state |0) and two e' g states |+) and |— 



where the local coordinate system of a perfect VC>6 octahedron was adapted [see Fig. DD(c)]. 
Note that the higher e g states d x 2_ y 2 and d z 2 remain degenerate. The recent polarized 
neutron diffraction study^ showed that, at each V site of a given V4 tetrahedron, only the 
dig level is occupied by an electron. 

The band structure and density of states calculated for the ferromagnetic (FM) state of 
Y2V2O7 by the DFT+U method are shown in Fig. [2](a) and Fig. EJ^a), respectively. Both the 
valence band maximum (VBM) and conduction band minimum (CBM) are majority-spin 
states. There is an indirect band gap of about 0.33 eV between the VBM state near the W 
point and the CBM state at the V point. This band gap is consistent with the experimentally 
measured activation energy 0.2 eV-. Our calculations show that the FM state of Y2V2O7 
has a nonzero band gap when U — J > 1.5 eV. The four states between —0.7 eV and eV 
(with zero set at the VBM level) are the occupied 3d states of the four V 4+ ions per unit cell. 
The electron density associated with the four states, plotted in the inset of Fig. EJ^a), clearly 
shows that these occupied orbitals are the a± g states, in agreement with experiment-^. The 
analysis of the partial density of states provides further insight into the electronic structures 
of Y2V2O7. For the spin majority part, the orbital levels are consistent with the trigonal 
crystal field splitting. For the spin minority part, the empty |+) and |— ) states have lower 
energies than the unoccupied |0) state. This is so because the intra-orbital onsite Coulomb 
interaction U is much larger than the interorbital onsite Coulomb interaction (£/'), and the 
trigonal crystal field splitting between the a\ g and e' g states. The spin down a\ g state lies 
slightly lower in energy than the spin down e g states, and both states are delocalized with 
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■ = 1/ \/3{d xy + d yz + d xz ) 
) = -1/V3{d xy + e 2 ^d yz + e-^ 3 d x ) 
> = l/V3(d xy + e- 2 ™/ 3 d yz + e 2m / 3 d xz ), 
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strong hybridization between them. The electron configuration of the V 4+ ion is shown 
schematically in Fig. [3]^b). 

The symmetric spin exchange parameters were extracted by mapping the energy differ- 
ences between ordered collinear spin states obtained from DFT+U calculations onto the 
corresponding energy differences obtained from the quantum Heisenberg Hamiltonian for 
the spin- 1/2 system: H = J2i<j Jij&i • Sj- We consider all symmetric spin exchange interac- 
tions up to the third nearest neighboring pairs [see Fig. QJb)]; the NN exchange J\ within 
each V4 tetrahedron and the next NN exchanges J2, J3 and J4. To evaluate these four spin 
exchange parameters reliably, we considered 33 different ordered spin states and then deter- 
mined them by performing a linear least-square fitting analysis^. Our calculations show that 
J x = -7.09 meV, J 2 = -0.07 meV, J 3 = -0.31 meV, and J 4 = -0.28 meV, namely, J 2 , J 3 
and J4 are negligibly small compared with the NN FM exchange J\. Using the calculated J\, 
we estimate the Curie- Weiss temperature 9 = zS(S + l)Ji/3k B = 6 x |(| + 1) Jx/3/cb = 122 
K, which is close to the observed value 106 K for Lu2V20t^. The calculated spin exchange 
parameters show that the magnetic ground state is the FM state, in agreement with ex- 
periment. Interestingly, we find that J\ is always ferromagnetic for any reasonable U and 
J values (U < 8 and J > 0), which is not in support of the previous prediction- that 
antiferromagnetism is favorable when J < 0.7. 

We now discuss the microscopic origin of the ferromagnetism in R2V2O7 by comparing 
the electron hopping processes between adjacent spin sites 1 and 2 for cases when they have 
the FM and AFM arrangements. By symmetry, the overall total energy gain by the hopping 
is twice that of the majority-spin a\ g state of site 1 (a^). In the FM case, the a\^ g state 
can hop to the majority-spin e' g and e g states but not a lg of site 2 according to the Pauli 
exclusion principle. In the AFM case, the a\^ g can hop to all minority-spin states of site 
2. The energy gain from the virtual hopping is larger for the FM case than for the AFM 
case because the empty up-spin states of site 2 are closer in energy to the filled a lg state 

1 "t" 

for the FM case and the hopping between a Xg state and the minority-spin a\ g state of site 2 
is negligible due to our finding that the minority-spin a\ g state lies much higher in energy 
than the minority-spin |+) and |— ) and a rather small transfer integral between a\ g states. 

Fig. [2]^b) shows the band structure of an AFM state in which there are two up- and 
two down-spins in every V4 tetrahedron [see inset of Fig. EJ^b)]. As in the FM case, there 
is an indirect band gap between the VBM state near the W point and the CBM state at 
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the T point. However, there is an important difference: The band gap of the AFM state 
is about 0.62 eV, which is almost twice that of the FM state. By taking the level of the 
Y 4s semicore level as the reference, we find that the VBM level of the FM state is at 
almost the same (only about 0.01 eV higher) as that of the AFM state. Therefore, the CBM 
of the AFM state is much higher than that of the FM state. The reason why the AFM 
state has a higher CBM and thus a larger band gap is illustrated in Fig. [3](d) . When the 
spins of two neighboring V ions have an FM coupling, the up-spin e' g states of the adjacent 
sites have the same energies so that the lowest energy state is lower than the e' g level by t', 
where t' is the hopping integral between the adjacent e' g states. In the AFM case, however, 
the lowest energy state is lower than the majority-spin e' g level by t' 2 /A, where A is the 
exchange splitting and A > if . Therefore, the AFM state has a higher CBM and thus a 
larger band gap. This fact naturally explains the negative MR observed just above the 
ferromagnetic Curie temperature T c . When the temperature is lowered towards Tq, the 
spins have a tendency to order ferromagnetically, but are not fully aligned. The application 
of an external magnetic field (about 5 T) will help align the spins ferromagnetically. Thus, 
the band gap decreases with increasing the magnetization, so that the resistivity of R2V2O7 
would be reduced by external magnetic field hence leading to the negative MR effect. This 
explanation is consistent with the observation that the maximum MR effect occurs at 75 K, 
right above the Curie temperature [See Fig. 3(a) of Ref.-]. We note that this new mechanism 
of negative MR should be also applicable to other ferromagnetic insulators (e.g., EuO). 

Let us now examine the magnetic properties that require the consideration of SOC. 
For two interacting spins, the SOC can induce DM antisymmetric interactions Hom = 
Djj • (Si x Sj) (Djj = Ddij, D = |Djj| and dy is a unit vector). According to the crystal 
symmetry, the DM vector for a V-V edge of each V4 tetrahedron is perpendicular to the 
V-V bond and is parallel to the opposite edge of the V 4 tetrahedron [Fig. Hfd)]. To evaluate 
the magnitude of the DM interaction term D, we consider two spin configurations shown 
in Fig. [He). In one spin configuration, the four spins are aligned along the X, Y, Z, and 
Z axes, respectively [see Fig. [TJ^a) for the definition of the global coordinate system XYZ\. 
From this configuration, we generate the other spin configuration by performing a reflection 
operation of each spin with respect to the XZ plane containing the spin site. The only 
difference from the first configuration is that the spin at the second site now points along 
the — Y direction. It can be easily shown that the two spin configurations have the same spin 
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exchange interactions and the same single-ion anisotropy interactions (see below). In terms 
of the quantum expression of the DM interaction Hamiltonian, it is found that the total 
DM interaction vanishes for the first spin configuration, but is Dy2/2 per V4 tetrahedron 
for the second spin configuration. By using the energy difference between the two spin 
configurations obtained from the DFT+U+SOC calculations, D is estimated to be 0.34 
meV. The D value is rather insensitive to the calculation parameters (U and J). Thus, D 
is of the order of 5% of the NN spin exchange J\ = —7.09 meV, i.e., \D/Ji\ = 0.048, which 
is almost an order of magnitude smaller than the value D/Ji = 0.32 deduced by Onose et 
al. from analyzing the magnetic field dependence of the thermal Hall conductivity in terms 
of their model for the magnon Hall effect. This raises a serious question as to whether the 
observed thermal Hall effect can be described solely in terms of the DM interaction. 

Another important consequence of the SOC interaction is that the magnetic moment 
of each spin site gets a preferred orientation in space with respect to the crystal lattice. 
However, it is commonly believed that a spin- 1/2 ion has no single-ion anisotropy because 
the spin doublet state is not split by the zero-field splitting term CSf = CI/4 for S — 1/2 (C 
is a constant and / is a 2 x 2 unit matrix) 1 ^. Here, we find that the spin- 1/2 V ions of R2V2O7 
have substantial single-ion anisotropy, which is not described by the usual zero-field splitting 
term. The single- ion anisotropy of the V 4+ ion can be estimated in two ways. First, we 
replace three of the four V 4+ ions in each V4 tetrahedron of Y2V2O7 with non-magnetic Ti 4+ 
Ti 4+ ions to obtain Y2TisVi07, which has no NN pairs of V 4+ ions. Our DFT+U+SOC 
calculations show that the easy axis of the spin- 1/2 V ion is along the three-fold rotational 
axis \z! in Fig. [1(c)] of the distorted V06 octahedron. The state with spin moments parallel 
to the z' axis is more stable than that with the spin moments perpendicular to the z'-axis 
by about 0.81 meV per V. Second, we consider two spin configurations to obtain a more 
accurate value of the single-ion anisotropy. In one configuration, all spins are along the 
easy axis directions such that, in a V4 tetrahedron, one spin is pointed out from the center 
and the remaining three spins are pointed to the center. From this spin configuration, we 
obtain the other spin configuration by rotating the directions of all the spins around the 
global Z axis by 120°. In the resulting spin configuration all spins are perpendicular to the 
easy axis directions. The two spin configurations are the same in the symmetric exchange 
interactions and in the DM interactions. Our DFT+U+SOC calculations for the two spin 
configurations show that the spin orientation along the easy axis is more stable than that 
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along the hard axis by 0.91 meV per V. Therefore, it is unequivocal that a spin-1/2 ion can 
have significant single-ion anisotropy, contrary to the general belief. It should be noted that 
single-ion anisotropy for a spin-1/2 ion is not excluded according to the general expression 
of the L • S Hamiltonian^ 1 ^. In the traditional effective spin Hamiltonian approach^ 1 ^, 
the energy for an orbital/spin basis function \LMl)\SMs) does not depend on the spin 
projection \Ms), whereas an exchange splitting is always present in any realistic magnetic 
system, as found in our first principles calculations. 

As described above, the single- ion anisotropy energy is much greater than the DM inter- 
action parameter in R 2 V207. We now estimate how the DM interaction parameter is affected 
when the single-ion anisotropy energy is neglected. The single-ion anisotropy for the spin- 
1/2 V 4+ ion cannot be described by any quantum spin Hamiltonian. Our calculations show 
that the single-ion anisotropy energy can nevertheless be written as H an i = AJ2i{&i • Z Q 2 
(|S| = 1/2 and A = —3.64 meV) when the spins are treated as classical vectors. Then, the 
total classical spin Hamiltonian for R2V2O7 is expressed as: 

H = ^2JSi-Sj + D^2 dy • (S, x Sj) + A ]T(S, ■ z'f. (2) 

i<j i<j i 

Using the spin Hamiltonian, our Monte Carlo simulation shows that the spin ground state 
is nearly FM with the moments aligned along the [001] direction, which is in agreement 
with the slight anisotropy displayed in the magnetization curved. If we renormalize the 
single ion anisotropy term into the DM interaction term, the Hamiltonian is rewritten as: 
H = J2i<j JSi " + D e ff J2i<j dij ■ (Si x Sj), where D e ^ is the effective DM interaction 
parameter to be obtained by neglecting the single-ion anisotropy. For two spin configurations 
that are related to each other by a mirror-plane symmetry [e.g., see Fig. QJe)], the energy 
differene can be used to extract D e ^ . We generate several random spin configurations to 
find that ^=p- can be as large as 20. This shows that the effective DM interaction can be 
much larger than the actual DM interaction when the single-ion anisotropy is neglected. 

In summary, R2V2O7 exhibits negative MR because its band gap depends on the spin 
arrangement, with the smallest gap for the FM state. The V 4+ ions of R2V2O7 exhibits 
an easy-axis single- ion anisotropy that is much stronger than the DM interaction term D, 
despite the common belief that spin-1/2 have no single-ion anisotropy. The effective D 
value evaluated can be unreasonably large when this anisotropy is neglected, as found from 
the analysis of the magnon quantum Hall effect of Lu 2 V 2 7 . Thus, the consideration of 
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the single-ion anisotropy is necessary to formulate a more complete theory for the observed 
magnon Hall effect. 
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FIG. 1: (a) The crystal structure of R2V2O7. The global coordinate system XY Z is indicated, (b) 
The spin exchange paths between V 4+ ions, (c) The two local coordinate systems used for the ideal 
VO"6 octahedron. The xyz coordinate system is defined for the ideal VOq octahedron with the x, y 
and z axes taken along the V-0 bonds. In the x'y'z' coordinate system, the z' axis is taken along 
one three- fold rotational axis of the ideal V06 octahedron. In R2V2O7 each VC>6 octahedron is 
axially compressed slightly with only one three-fold rotational axis, (d) The DM vectors Djj with 
i < j of the V tetrahedron, where i and j denote the V site labels, (e) The two spin configurations 
used to extract the DM parameter, (f) The two spin configurations used to extract the single-ion 
anisotropy parameter. 
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FIG. 2: (color online) (a) The band structure of the FM state from the DFT+U calculation, (b) 
The band structure of the AFM state (with two up- and two down-spins in each V4 tetrahedron) 
from the DFT+U calculation. Solid lines and dashed lines represent up-spin and down-spin bands, 
respectively. 




FIG. 3: (color online) (a) The total and partial density of states of the FM state from the DFT+U 
calculation. The inset shows the density distribution associated with the occupied V 3d 1 states, 
(b) A schematic illustration of electron configuration of the V 4+ ion. (c) The interactions between 
the 3d-states of two adjacent V 4+ ions when their spins have the FM and AFM arrangements. 
Only the hopping processes which lead to the overall energy lowering of the up-spin a\ g state of 
the spin site 1 are shown, (d) The interactions between the 3d-states of two adjacent V 4+ ions 
leading to the CBM positions when their spins have the FM and AFM arrangements. 
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